Communication/ in Information Science and management engineering 



cisme 



Simulation and Implementation of BPSK BPTC and 
MSK Golay Code in DSP Chips 

Yihua Chen, Juehsuan Hsiao, PangFu Liu, Kunfeng Lin 

Oriental Institute of Technology Institute of Information and Communication Engineering New Taipei city, Taiwan 
ff011@mail.oit.edu.tw, ff014@mail.oit.edu.tw, mars7699@gmail.com 



Abstract-This study discussed the decoding mechanism of block 
product turbo code (BPTC) (196,96) and simulated its efficiency, 
used the Hamming (15,11) and Hamming (13,9) block channel 
code combinations and block interleaving to construct a BPSK 
modulation and BPTC (196,96) coding system in the concept of 
feedback encoding in turbo code. The decoding mechanism was 
meticulously described to explain the advantage of this turbo 
feedback structure. This study also analyzed and simulated 
special combination of MSK modulation with Golay (20, 8) 
channel code. Golay (20,8) decoding correction capability was 
comprehensively discussed. The simulation results indicate MSK 
Golay (20,8) is less robust than BPSK Golay (20,8) for 2.8dB 
against noise, but BPSK BPTC (196, 96) is 3dB more robust than 
MSK Golay (20, 8) and it is also a better choice than BPSK Golay 
(20, 8), if signal Ej/N is stronger than 5dB. Whenever the C code 
system program is constructed, the system program should be 
optimized and realized on the DSP platform. Both MSK Golay 
and BPSK BPTC C program optimization efficiency statistics are 
analyzed. The implementation can be applied to the wireless 
communication system of SDR concept, so that different 
modulation and encoding modes can adjust parameters in 
firmware directly according to different system specification 
requirements. 

Keywords-component; BPSK; MSK; Golay Code; BPTC; 
Coding Gain; Syndrome; Perfect Code; Error Pattern; BER; SDR; 

I. INTRODUCTION AND BACKGROUND 

When a wireless communication system is in a non-Line- 
of-sight (NLOS) environment, the signal will be refracted, 
reflected, diffracted and scattered when passing through 
different features in the course of transmission. For a real-time 
communication system, how to reduce the error probability of 
digital signal in propagation is a very important issue. The 
channel coding is a technology that reduces the signal 
distortion resulted from multi-path interference and AWGN 
noise in receiving. Since the channel coding must add excess 
bits to the original data in transmitting (except Trellis Coded 
Modulation, TCM), called parity bit, so that we can detect and 
correct the error bits of signal. The additional parity bits will 
reduce the energy of received signal symbol, so that the symbol 
error rate in receiver increases. However, if more error bits can 
be corrected to compensate for the symbol error rate increased 
by the decrease in symbol energy, the bit error rate of the 
whole received signal still can be reduced. If only the same bit 
error rate is required to be kept, the power efficiency Eb/NO of 
the received digital signal after decoding can be reduced. In 
other words, the coding gain is increased [1], the precondition 
is that the bit error rate must be identical. 

G(dB) = (E t /No)uncoding(dB) - (E b /N ) coding (dB) 



To attain the coding gain, the signal-to-noise (SNR) of the 
received signal must reach a certain Eb/NO threshold, 
otherwise the addition of coding parity bits increase the symbol 
error rate of the received signal. The decoding mechanism 
cannot compensate for this loss, and there will be no coding 
gain at all. Therefore, the selection of coding rate is very 
important; in other words, adopt different encoding modes and 
different coding rates should be adopted according to different 
system specification requirements. The coding rate is 
determined by the number of information bits and parity bits 
that influences the decoding capability of channel code 
significantly. 

The bit error rate of signal in transmission can actually be 
improved by adding channel coding in a wireless channel 
environment. The excess bandwidth is taken out from limited 
bandwidth for channel coding to make the receiver have lower 
bit error rate at the same transmitting power. Alternatively, the 
required communication quality can be attained at lower 
transmitted power, from which we have coding gain of the 
channel code. Both benefits can only be achieved when the 
Eb/NO threshold of received signal are reached. The Eb/NO 
threshold effect will be discussed in this study. 

In a wireless channel environment, the most frequent 
problem is the interference and noises from different 
environments, such as atrocious weathers, obstacles and user 
moving. Therefore, continuous burst errors are likely to happen, 
that will cause great difficulties in coding. Therefore, we must 
use block interleaving to break up the burst errors into 
discontinuous random errors, so as to exert the capability of 
channel coding to correct errors. 

Generally, adding more redundancy bits can correct more 
random error bits, but the difficulty in coding and decoding 
will increase, and the complex decoding algorithm requires a 
larger memory and complex software and hardware design [2]. 

The turbo code was introduced by Berrou, Glavieux and 
Thitimajshima in 1993, published in literatures [3] [4] which 
indicate that this scheme can reach the bit error probability 
of 10-5 in additive white Gaussian noise channel and Eb/No of 
0.7 dB BPSK modulation at 1/2 code rate. This encoding mode 
combines more than two simple block codes or convolutional 
codes for larger coding gain. The concept behind turbo 
decoding is to transmit the soft decision from the output end of 
the first decoder to the input end of another decoder, and iterate 
this process for several times, so as to generate relatively 
reliable decisions and add in interleaving coding. 

In 1948, Hamming put forward an important error- 
correcting code, named Hamming code [5]. It uses parity check 
matrix to detect and correct errors; however, its detection and 
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correction capability is limited, it can only detect 2-bit errors or 
correct 1-bit errors. 

This paper discusses the block product turbo code (BPTC) 
which is classified as one of block turbo code concatenation 
forms. The Hamming code can detect two-bit error or correct 
one bit error. The BPTC uses two Hamming codes for 
"column" coding and "row" coding, it has improved the 
Hamming code correcting only one error. In addition, the 
BPTC carries out block interleaving coding for disorganizing 
the transmission sequence before transmission, so as to avoid 
burst errors when the signal meets multi-path channel in the 
channel. This paper will also explain the principle of error bit 
correcting mechanism of BPTC channel coding. 

There are six sections in this paper. The second section 
introduces the operating principle of the two Hamming codes, 
Hamming (15, 11) and Hamming (13, 9), used in BPTC code. 
The third and fourth sections discuss the BPTC and Golay (20, 
8) code design principles. In the fifth section, we discuss all the 
simulated data of BPSK BPTC and MSK Golay. In the last 
section we infer some conclusions of these results. 

II. HAMMING CODE 
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Among which, Vi~Vk are linearly independent vectors that 
can generate all code vectors. The data of the transmitting 
terminal are usually expressed in column vectors, therefore, the 
sequence of k message bits, i.e. the message m is expressed as 
l*k matrix. If we use matrix notation to indicate the generation 
process of code word, i.e. the product of m and G, 

U = mG 

The two Hamming codes used in this paper are Hamming 
(15, 11) and Hamming (13, 9), as defined above, the (15, 11) 
code means n is 15 and k is 11, the redundancy bit added in the 
coding process is 4, and the correction capability is 1. The (13, 
9) code means n is 13, k is 9, the redundancy bit added in the 
coding process is 4, and the correction capability is 1. 

The generator matrices are shown in Fig. 2 and Fig. 3. 



A. Hamming Encoding 

The Hamming code is a sort of binary linear system block 
code (as shown in Fig. 1). Hamming put forward an important 
single-error-correcting code, using parity check matrix (H) to 
detect and correct errors. Hamming code is an important 
forward error correction (FEC) in theory and practice so far. 
Furthermore, Hamming code (7, 4) and (16, 11) is very limited 
perfect code in systematic code. But is not for Hamming (15, 
11) and Hamming (13, 9), used in this paper. 
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Figure 2. Generator matrices Hamming (15, 11) 



Figure 1. Schematic diagram of Hamming coding process 

Hamming code is a simple type of systematic code, 
described as the following structure. 

Blocklength:n = 2 m -l 

Number of data bits: k = 2 m -m-l 

Number of check bits: n-k = m 

Minimum distance: d min = 3 -> Correct single bit error 

(n,k) = (2 m -l,2 m -l-m) 

Generally speaking, we can use the following k x n array to 
define the generator matrix (G) 
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Figure 3. Generator matrices Hamming (13, 9) 

B. Hamming Decoding 

In order to decode the received signals, we need to define a 
parity check matrix and a syndrome. There is a (n-k)xn matrix 
H in each generator matrix G, so that the columns of G are 
orthogonal to the columns of H, i.e. GH T =0, the H T is the 
transpose matrix of H. In order to meet the orthogonality of 
system coding, the component of matrix H can be expressed as 



H = [I n _! 



P T ] 



Therefore, the matrix H can be expressed as 
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H T = 



U is the code word derived from matrix G if and only if 
UH T = 0. Let r be the vector received by the receiving terminal, 
so the r can be described as 

r = U + e 

Among which, e is the error vector resulted from the 
channel. There are 2 n -l error patterns in nonzero position in the 
space formed by 2 n value groups. The syndrome is defined as 

S = rH T 

The syndrome is the result of parity check implemented in r, 
judge whether r is an effective element in the codeword set. 
Based on development of equation 



S = (U + e) H T 



UH T + eH T 



However, for all elements in codeword set UH T =0, 
therefore 

S = eH T 

Since the correction capability of Hamming(15,ll) and 
Hamming(13,9) is 1, meaning the error pattern is one selected 
from n, so the syndrome tables can be set up, the syndrome 
tables are shown in Fig. 4. 
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(a.) Hamming (15, 11) 



(b.) Hamming (13, 9) 



Figure 4. Hamming code syndrome tables 

The error correction coding procedure is as follows: 

1. Use S = rH T to calculate the syndrome of r 

2. Find out common first item (error pattern) e j} its 
syndrome equals rH T 

3. This error pattern is supposed to be the error caused 
by channel 

4. The identified receive correction vector or code word 
equals U = r + e 7 



III. BLOCK PRODUCT TURBO CODE design principle 

The 96 messages must be arranged in 11x9 matrix form 
before BPTC coding, the 3 absent messages are set as 0, placed 
in the last three positions of the matrix (as shown in Fig. 5), 



Step 3. Hamming(15,11) encoder 








Step 1 . Message [96] 


BRUM 


BMC) 


H.KO.'L 


B_R<H0l 


l(IC) 


«»t 


hsi 


IO 


CSf 


W*t 


WJ 


WJ) 


1C> 


1(1) 


1(5) 


H_RJ|1| 


H_RI,1, 


H_R1|1| 


H_R1!»| 


l(=l| 


WB 


1(1*) 


1") 


«H» 


I(Uf] 


Hi*] 


1(11) 


nam 


«3JJ 


WD 


H_W|i| 


H_Ki(!) 


11 R'.'l:. 


H_RS((H 


K«t 


ran 


H«> 


Ul*> 


K5») 


i :-, 


K!*) 


K»> 


!(!*> 


K!J) 


k::> 


H.KW 


H BJi:, 


B.RHI) 


Bjam 


««» 


[(«j 


U-U) 


K«i 


CIS) 


WM1 


US -J 


HJSJ 


Ipi) 


1(J4> 


1(31) 


KRMSl 


H_R*iii 


K_l»|l| 


H_RJ(»| 


l(«| 


p«j 


K«) 


n*i! 


U») 


K'-M) 


H43i 


1(J7> 


]('*> 


1(JS> 


KU) 


n_*m 


H KS,:i 


H.K^It 


Bjtm 


K*ft 


W*) 


K«) 


!,,:, 


l!»l> 


1 ,.,.. 


KW) 


KM> 


CiS- 


It!*) 


Kff) 


H_R£|1» 


H_M(i) 


ttjtmt 


H R( 


1(7*1 


PTJ) 


K7-U 


K7i| 


(paj 


IPIJ 


a :o) 


W*»> 


]£«» 


1(47| 


I(Mt 


H.R^it 


H_m;i 


H.RT(I» 


H.R-(0| 


na- 


HS6) 


].Sf. 


K*J> 


nsi) 


LIS!) 


HSli 


USD, 


l£T») 


1[7*> 


H"> 


H Rft|J| 




H s.i . 


il tam 


na 


K.t 


K(4| 


l(M> 


K*u 


■ONI 


H*:> 


Il«l 


1C**> 


1(«) 


[,v, 


H_C1^«) 


H_CIJ«H 


h_ci;c«) 


H_t U|t> 


E_Cl*l«) 


E_CSH«) 


H_Ct((i) 


e_c im 


E_C*(ti) 


a_cm 


H_<JiOi 


H_CJ(»h 


H_C3(*( 


H_Cl(t) 


H_t0(»h 


H_t]«l] 


H_CIHD 


h < i:i. 


H i 11. 1 


H t I'lli 


H_C«i) 


H_C5(l) 


H_C-T(1) 


E_C*(l) 


H *?t1i 


K * ul, 


H_CJ(1> 


B_c:(i> 


B_C1(1) 


H CO. 11 


HJtMZ) 


h.cim:) 


h ! i: :, 


B_CLl(i) 


h ( !■■,:, 


B_CKi) 


B.CSCI) 


B.C7C) 


E_«C) 


a.cfo 


H_C«!} 


H.CKi) 


B.C!C) 


B.C1C) 


h ci.,:. 


H_C1-H!) 


H_tlJ(J» 


H_ti:iJ) 


H_t U(l| 


H_Cl»U:i 


H_CJHJ) 


H_CS(S) 


H_C7<Jj 


E C*(J) 














Step 2. Zero padding [3] 








Step 4. Hamming(13,9) encoder 



Figure 5. BPTC (196, 96) Matrix 

which are R(0), R(l) and R(2), one R(3) is added in the coding 
process for replacing the position during interleaving after 
coding, so it is not placed in the messages to be encoded. 

A. BPSK BPTC Communication System Architecture 

The digital signal will be coded by BPTC (196, 96) and 
modulated by binary phase shift key (BPSK) modulator before 
transmission. When the signal passed through the AWGN 
channel, the received signal is demodulated by BPSK 
demodulator and sent to the decoder of BPTC (196, 96). The 
system block diagram is shown in Fig. 6. 
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Figure 6. BPSK BPTC system block diagram 

B. BPTC Encoding 

The Hamming (15, 11) is encoded according to the above 
Hamming encoder, I (0) to I (10) are a set of messages, H_R0 
(0), H_R0 (1), H_R0 (2), H_R0 (3), I (11) to I (21) derived 
from encoding are a set of messages, H_R1 (0), H_R1 (1), 
H_R1 (2) and H_R1 (3) are derived from encoding, till I (88) 
to R (2) are a set of messages, afterwards, Hamming (15, 11) 
encoding is finished. 



CISME Vol.1 No.4 2011 PP.46-54 www.jcisme.org ©World Academic Publishing 

-48- 



Communication/ in Information Science and Management Engineering 



cisme 



After Hamming (15, 11) was encoded, the message array 
has changed to 15x9 form, and then Hamming (13, 9) will be 
encoded. 

Encode Hamming (13, 9), I (0), I (11), I ( 22),..., I (88) are 
a set of messages, H_C0 (0), H_C0 (1), H_C0 (2), H_C0 (3) 
are derived from encoding, I (1), I (12), I (23),..., I (89) are a 
set of messages, H_C1 (0), H_C1 (1), H_C1 (2) and H_C1 (3) 
are derived from encoding, the rest can be deduced accordingly, 
till H_R0 (3), H_R1 (3), H_R2 (3),..., H_R8 (3) are a set of 
messages, afterwards, Hamming (13, 9) encoding is finished. 

When the row and column encoding is completed, the 
matrix has changed to 15x13 form, and then the interleaving is 
carried out, due to the interleaving needs the sequence number, 
the encoded matrix must be renumbered, H_R0 (0) should be 
the 12th message, H_R0 (1) is the 13th message, and I (11) will 
be the 16th message, the rest can be deduced accordingly, the 
results are shown in Fig. 7. 
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Figure 7. Interleaving indices for BPTC (196, 96) 

The index in Fig. 7 is the redefined message number, 
multiply the message number by 13, divided by 196, take the 
remainder, the new interleave index number can be formulated 
as the following statement 

Interleave Index = Index x 13 modulo 196 



The result of interleaving can be obtained from the new 
numbers derived from the above equation; the messages should 
be transmitted according to the new number sequence. 

It should be paid special attention to that the new number of 
R (3) which was not placed in the messages originally is placed 
in the position of 183 after interleaving, and the I (11) in the 
original messages is placed in the position of 195 after 
interleaving, only the messages of new numbers from to 194 
will be sent out, the position of 195, i.e. the original I (11) will 
not be transmitted together. 

C. BPTC Decoding 

The messages received by the receiving terminal are in the 
form of 15x13 matrix, de-interleaving must be carried out 
before decoding to restore the original message sequence. 
Decoding can be carried out after de-interleaving, the decoding 
should be in such an order that encoded late should be decoded 
first, that encoded early should be decoded late. 

Decode Hamming (13, 9) according to the aforesaid 
Hamming decode, I (0), I (11), I ( 22),..., H_C0 (2) and H_C0 
(3) after sequence restoration are a set of messages, decode and 
correct I (0), I (11), I (22),..., I (88), I (1), I (12), I (23),..., 
H_C1 (2) and H_C1 (3) are a set of messages, decode and 
correct I (1), I (12), I (23),..., I (89), the rest can be deduced 
accordingly, till H_R0 (3), H_R1 (3), H_R2 (3),..., H_C14 (2) 
and H_C14 (3) are a set of messages, afterwards, Hamming (13, 
9) decoding is finished. 

After Hamming (13, 9) decoding, the message array has 
changed to 15x9 form, and then Hamming (15, 11) will be 
decoded. 

Decode Hamming (15, 11), I (0), I (1), I (2),..., H_R0 (2) 
and H_R0 (3) are a set of messages, decode and correct I (0), I 
(1), ..., I (10). I (11), I (12), I (13),..., H_R1 (2) and H_R1 (3) 
are a set of messages, decode and correct I (11), I (12), ..., I 
(21), the rest can be deduced accordingly, till I (88) to H_R8 (3) 
are a set of messages, afterwards, Hamming (15, 11) decoding 
is finished. 

D. BPTC error correcting mechanism analysis 

Although the BPTC is composed of quadratic Hamming 
code, its probability of being corrected in the second dimension 
coding is increased by using the fundamental characteristics of 
turbo code. 

As shown in Fig. 8 (a), if the column vectors only have 
error in position (1.), they can be certainly corrected by using 
the characteristics of Hamming. When the column vectors have 
error in position (2.), if the row vectors of the second 
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dimension have no error, the error in position (2.) can be 
corrected as well, the key is whether the second dimension has 
error in position (3.), the same situation can be observed in any 
column, therefore, return to the column vectors, although the 
first dimension can correct only one message, the other 
messages can be corrected in the second dimension probably. 
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Figure 8. Schematic diagram of two-dimension correction capability 

For example, in Fig. 8 (b), when decoding the column 
vectors of the first dimension, positions (1.), (2.), (3.), and (4.) 
have only one error message, so these errors will certainly be 
corrected to the original messages. When decoding the row 
vectors of the second dimension, due to positions (5.) and (6.) 
can be probably corrected in the first dimension decoding, 
there will be three situations at this point, the first one is when 
the position (5.) is corrected in the first dimension decoding, 
the position (6.) will be corrected in the second dimension 
decoding. The second situation is when the position (6.) is 
corrected in the first dimension; the position (5.) will be 
corrected in the second dimension. The third situation is when 
the positions (5.) and (6.) cannot be corrected in the first 
dimension, the row vector decoding of the second dimension 
will still correct the errors in positions (5.) and (6.), for row 
vectors, (5.) and (6.) are one single error message, they can be 
corrected certainty. 

E. BPSK BPTC C program optimization efficiency analysis 
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Figure 9. BPSK BPTC optimization efficiency statistics 

After calculation, the fundamental operations and 
functional equations are shown in Fig. 9., the double represents 
floating-point number operation, » and « represent bit right 
shift and left shift respectively, compare instead of judgment 
operation and the compare represents identification operation. 
Most of quadratic square operations are changed into right shift 
and left shift, due to one left shift of number 1 in C code is 
equivalent to first power of 2, but the efficiency of left shift and 
right shift is much better than that of power, such conversion 
can improve the computing speed. 



IV. GOLAY (20,8) CODE DESIGN PRINCIPLE 

This paper also discusses Golay (20, 8) which is one type of 
block code. The minimum Hamming distance of general Golay 
(23, 12) code [6] [7] d min = 7, number of corrected error bits t = 
3. All the 3-bit errors can be corrected, so that Golay (23, 12) 
code is the so-called perfect code [8], for it meets the Hamming 
bound of the following equation 



r - k= £ 



' r n \ 



The extended non-perfect Golay (24,12) code emphasizes 
the relation among n, k and t, namely, the coding rate can have 
simple integer ratio for convenient digital circuit 
implementation. 

For the Golay (20, 8) code discussed in this paper, code 
word length after coding n = 20, the number of information bits 
k = 8, it is of extended Golay code. The minimum Hamming 
distance of this Golay (20, 8) code is increased to d min = 8, 
besides correcting all 3-bit errors, the error correcting 
capability is improved to partial 4-bit errors. However, even the 
larger minimum Hamming distance d min = 8 is not emphasized, 
the correction capability is improved by increasing the 
minimum Hamming distance, so as to reduce the signal error 
probability. 

The predominance of Golay (20, 8) code is the large 
corresponding syndrome list, therefore, more error patterns can 
be corrected, certainly its decoding circuit is complicated 
relatively. The operating principle and simulation result of 
Golay (20, 8) code will be detailed below. 

In the following equation, g(x) is the generator polynomial 
of Golay (20, 8) code. 



9(x) 



: x 11 + x 10 + x 6 + x 5 + x 4 + x 2 + 1 



6 165 8 



The generator matrix in Fig. 10 can be divided into two 
parts, the left part is an 8^8 identity matrix, and the right part is 
an 8x12 parity matrix. 
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Figure 10. The generator matrix of Golay (20, 8) 

A. Golay (20, 8) encoding 

The information of the original data of Golay (20, 8) is 8- 
bit, while added with 12-bit parity bit for coding; it changes 
into 20-bit code word transmission after coding. The 
transmitting terminal uses generator matrix for coding, the 1x8 
information data matrix is multiplied by 8x20 generator matrix 
in Table 1 to obtain 1x20 coding code word matrix. 
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B. Golay (20, 8) decoding 

The receiving terminal uses parity check matrix for 
decoding, the structure can be set as 



H = [/„_„ I P T ] 



Thus, the columns of generator matrix (G) are orthogonal to 
the columns of parity check matrix (H), i.e. 

GH T =0 

Therefore, HT can be expressed as 



H T = 



I 



Next, the decoding is divided into four steps: 



Stepl. Multiply received code word (1x20) by H T to obtain 
(1x12) syndrome vector, 

Step2. The only error pattern was educed out by the 
corresponded syndrome. The vector matrix of error pattern is 
(1x20), 

Step3. Add this error pattern to the received codeword, so 
as to obtain the corrected code word, 

Step4. Since the Golay (20, 8) code is of systematic code, 
we can take out the original information from the corrected 
code word, see the following equation, 

[received codeword] i x2 o = [Mi M 2 ...M 8 Pi P 2 P 3 ...P12 ] 

1x20 

therefore, the decoded information is [Mi M 2 M 3 . . .M 8 ]. 

C. MSK Golay communication system architecture 

The digital signal will be coded by Golay (20, 8) and 
modulated by minimum shift key (MSK) modulator before 
transmission. When the signal passed through the AWGN 
channel, the received signal is demodulated by MSK 
demodulator and sent to the decoder of Golay (20, 8). The 
system block diagram is shown in Fig. 11. 
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coding (here also announced as int type), it is a 1x20 block 
codeword after coding, and then enters MSK modulator (here 
the data type is int input, double output), AWGN channel (data 
type is equivalent to MSK modulator), MSK demodulator (data 
type is equivalent to MSK modulator), the received signal of a 
block information bits of 1x8 bits is changed into an int output 
after decoding. 

The bit error rate (BER) is figured out by comparing the 
original transmission signal with the finally decoded received 
signal. 

The program must be optimized in order to reduce 
operation steps and simplify the functional equation. The Box- 
Muller [9] is used for the generation principle of AWGN noise, 
the equation is as follows 



A /-21n(L/ 1 )cos(2^-L7 2 ) 

Among which, Ui and U 2 are random number of to 1, we 
should use the program produced cos function look-up table 
rather than the C language built-in function, so as to increase 
the computing speed for program optimization. 

This paper uses a simpler equation to replace cos function 
for program optimization. The equation is as follows. 

2.0*rand()/RAND_MAX-1.0 

Among which, rand()/RAND_MAX is a random number of 
to 1, RAND_MAX=32767, image this range to the range of 
cos function, changes the original range of to 1 to -1 to 1. 

Fig. 12 is AWGN in Rayleigh probability density function 
(PDF) profile, the blue solid curve is theoretical value, the red 
dotted curve is the simulation result of C language optimized 
program, here, the SNR is 20dB. 



0.012 




Figure 12. The AWGN in Rayleigh probability density function profile 

D. Golay (20,8) encoding processes 

The following (l x 8) matrix is the data information. 



[IX 



10 10 



Figure 11. MSK Golay system block diagram 



Multiply a set of data information of matrix (l x 8) by 
generator matrix in Fig. 9 to obtain the coded code word (1x20) 



First, in C-code program the transmission signal of a block 
of 1x8 bits is changed into an int, enters the Golay encoder for 
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E. Golay (20,8) decoding processes 

The following 20x12 transposed parity check matrix (H T ) is 
used for decoding. 
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Stepl. Multiply the coded code word by transposed parity 
check matrix to obtain syndrome (1x12) 



| 1 1 1 1 1 lj 

Step2. List the syndromes corresponding to error patterns 

From the look-up table of t = 1,2,3 [10], each of the 1350 
error patterns corresponds to a unique syndrome, in other 
words, 1-bit to 3-bit errors can be corrected completely. The 
original error pattern is that the number of combinations of four 
error bits in twenty code word bits is 4845. However, as 12 bits 
are a group in syndrome, there are 2 12 = 4096 syndromes. 
Therefore, when the information has 4 bit errors, only part of 
them can be corrected. 

Step3. Use the syndrome list to find out the corresponding 
error pattern, add the error pattern to code word, so as to 
correct 1-bit to 3-bit errors and partial 4-bit errors and reduce 
the error probability. 

F. MSK Golay C program optimization efficiency analysis 



respectively, the compare represents judgment operation of if, 
switch or ?. 

According to this table, power operation is avoided when 
composing the program, the power of 2 is replaced by left shift, 
for example 1 « 1 is 2 to the first power, 1 « 2 is 2 to the 
second power, thus, the program can be optimized and realized 
on DSP easily. 

V. COMPARISON OF C CODE ALGORITHM 
SIMULATION 

When we simulate a communication system based on 
BPSK modulation and channel AWGN [11], and keep the error 
probability at 10" 5 , the power consumed is 9.5 dB (as shown in 
Fig. 14). Hamming (15, 11) and BPTC (196, 96) are compared 
in the same AWGN channel environment, disregarding the 
excess bandwidth consumed by coding, if the error probability 
should be kept at 10" 5 as well, the BPTC has about 0.6 dB 
coding gain compared with Hamming code. 

If the Ej/Nq must be kept at 8 dB in the same AWGN 
channel environment, the error probability of BPTC is 
0.000006140625, and the error probability of Hamming (15, 11) 
is 0.0000299, it means the error probability of BPTC is lower 
than that of Hamming code by about 10" 1 at the same 
transmission power. 
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Figure 14. BPTC and Hamming error probability curves 

Fig 15 shows the MSK error probability curves calculated 
by theoretical MSK (blue solid) and simulated by C code (Red 
dotted). When P b = 10" 5 , the theoretical MSK E^Nq is at 
12.6dB, and the C code MSK Ej/No is at 12.5dB, obviously, the 
two curves are approximately identical. 
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Figure 13. MSK Golay optimization efficiency statistics 

Fig. 13 shows the fundamental operations and functional 
equations used in the system program after optimization, the 
double represents the floating-point number operation, » and 
« represent bit right shift and left shift respectively, log and 
sqrt represent log and square root in mathematical operation 
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Figure 15. MSK error probability curves simulated by C code (Red dotted) 
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Fig 16 shows the error probability curves without coding 
and with Golay (20, 8) coding, when P b = 10" 5 , the uncoded 
Ej/Nq is about 12.5dB, the Golay (20, 8) of t = 4 is about 
11.2dB. It is observed that at the same error probability, the 
coded Ej/No is decreased by 1.3dB, the coding gain is 1.3dB, so 
the advantage of taking Golay (20, 8) is quite obvious. 

Afterwards, the Golay (20, 8) of t = 4 in Figure 15 is 
compared with the Golay (20, 8) curve with correction 
capability up to t = 3 only. At P b = 10" 5 , Golay (20, 8) of t = 4 
is about 11.2dB, and the Golay (20, 8) of t = 3 is 12.2dB, the 
difference is ldB. Therefore, the Golay (20, 8) of t = 3 has not 
reached the due code efficiency. Because there should be 4096 
syndromes for correcting 4096 error patterns correspondingly, 
actually, there are only 1350 syndromes when t = 3, the error 
patterns that can be corrected are decreased by 67.04%. 
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Figure 16. The error probability curves without and with Golay (20, 8) code 

In addition, the Et/N thresholds of BPSK Golay (20, 8) of t 
= 4 (green and orange) and MSK Golay (20, 8) of t = 4 (pink 
and brown) are compared in Fig 17. The Ej/Nq threshold of 
BPSK Golay (20, 8) is at 5.2dB; nevertheless, MSK Golay (20, 
8) is at 8dB. There is a 2.8dB difference between them. That 
indicates BPSK Golay (20, 8) is more robust than MSK Golay 
(20, 8) for 2.8dB against noise. At same bit error probability Pj, 
= 10" 5 , MSK Golay (20,8) of t = 4 Ej/No is 11.2dB, yet, BPSK 
Golay (20, 8) of t = 4 Ej/No is 8.4dB. There is a same 2.8dB 
difference between them. 

With regard to Hamming (16, 11, 4) code, there is 1.2 dB 
lower Ei/No threshold compared to Golay (20, 8) of t = 4 in 
both BPSK (4dB) and MSK (6.8dB). According to this result, 
we infer that Hamming (16, 11, 4) is more solid than Golay (20, 
8) of t = 4 against noise and getting better performance in 
larger Efa/iVo. 

At last, we compare the E^No thresholds of BPSK BPTC 
(196, 96) (red and pink) and MSK Golay (20, 8) of t = 4 (dark 
blue and brown) in Fig 18. The Et/N threshold of BPSK BPTC 
(196, 96) is at 5dB; nevertheless, MSK Golay (20, 8) of t = 4 
is at 8dB. There is a 3dB difference between them. That 
indicates BPSK BPTC (196, 96) is 3dB more robust than MSK 
Golay (20, 8) of t = 4 and only 0.2dB more robust than BPSK 
Golay (20, 8) of t = 4. From Fig. 18, we can also observe that 
BPSK BPTC (196, 96) is getting better performance when the 
signal Eb/No exceeding the cross over point and approaching 
0.5dB gain relative to BPSK Golay (20, 8) of t = 4. 
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Figure 17. Comparison of the EM threshold of BPSK Golay (20, 8) of t = 4 
(green and orange) and MSK Golay (20, 8) of t = 4 (pink and brown) 
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Figure 18. Comparison of the EM threshold of BPSK BPTC (196, 96) (red 
and pink) and MSK Golay (20, 8) of t = 4 (dark blue and brown) 

VI. CONCLUSIONS AND FUTURE PROSPECTS 

This paper observes the efficiency of BPTC, and compares 
the complexity and coding gain with Hamming code, and the 
physical mechanism of decoding efficiency difference between 
Hamming code and BPTC. 

We can infer in this study that BPSK BPTC (196, 96) is a 
better choice than BPSK Golay (20, 8) of t = 4, if signal Ej/Nq 
is stronger than 5dB; nevertheless, BPSK Golay (20, 8) of t = 
4 will be a better choice than BPSK Hamming (16, 11, 4). 

In order to improve the communication quality at limited 
transmission power, it is observed from Figure 16 in this paper 
that when Ej/Nq is at 8dB, the uncoded MSK P b is 0.00519; 
MSK Golay (20, 8) is 0.00526. Once the transmission power 
reaches the certain threshold, the Golay (20, 8) t = 4 coded 
error probability curves have obvious difference with the 
uncoded one. For example, when a communication quality is 
required P b = 10" 5 , the uncoded MSK E^/Nq is 12.5dB, the 
MSK Golay (20, 8) t = 4 is 11.2dB. Therefore, the coded Et/N 
is decreased by 1.3dB and the coding gain is 1.3dB in the same 
communication quality. 

Although, the Et/N threshold of BPSK Golay (20, 8) is 
lower than the Ej/Nq threshold of MSK Golay (20, 8), from 
Figure 17, we can see significantly that both coding benefits of 
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MSK Golay (20, 8) of t = 4 and BPSK Golay (20, 8) of t = 4 
can only be achieved when the Ef/No threshold of received 
signal are reached. 

MSK modulation scheme has a great advantage of power 
efficiency in peak to average power ratio (PAPR), but less 
efficient in error performance compared to BPSK modulation 
scheme. At same bit error probability P b = 10" 5 , uncoded BPSK 
Ej/Nq is about 9.6dB, there is a approximate 2.9dB gain 
compared to uncoded MSK. At the same power efficiency, 
Ej/Nq = 8.4dB, BPSK Golay (20, 8) of t = 4 can reach P b = 10" 5 , 
but, MSK can reach worse performance than P b = 10~ 3 . 

Besides the aforesaid discussions, this paper uses 4096 
syndromes of Golay (20, 8) when d min = 8, compares correcting 
partial 4-bit errors and correcting only 3-bit errors, it is 
observed that the error patterns that can be corrected are 
increased by 67.04% when 2 12 syndromes are used completely. 

This paper uses C language for simulation analysis, when 
the C code system program is constructed, the system program 
is optimized and realized on the DSP platform, fast verification 
and SDR concept can be realized by replacing modules, such 
as modulation modules, channel modules and coding modules, 
so that different modulation and encoding modes can adjust 
parameters in software directly according to different system 
specification requirements, so as to economize hardware cost 
and product development time. 

The further goal of this implementation is to be applied to 
the wireless communication system of SDR concept, so that 
different modulation and encoding modes can adjust 



parameters in firmware directly according to different system 
specification requirements. 
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